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ABSTRACT 

The antizyme protein, Oaz1, regulates synthesis 
of the polyamines putrescine, spermidine and 
spermine by controlling stability of the polyamine 
biosynthetic enzyme, ornithine decarboxylase. 
Antizyme mRNA translation depends upon a 
polyamine-stimulated +1 ribosomal frameshift, 
forming a complex negative feedback system in 
which the translational frameshifting event may be 
viewed in engineering terms as a feedback control- 
ler for intracellular polyamine concentrations. In this 
article, we present the first systems level study of 
the characteristics of this feedback controller, using 
an integrated experimental and modeling approach. 
Quantitative analysis of mutant yeast strains 
in which polyamine synthesis and interconver- 
sion were blocked revealed marked variations in 
frameshift responses to the different polyamines. 
Putrescine and spermine, but not spermidine, 
showed evidence of co-operative stimulation of 
frameshifting and the existence of multiple ribo- 
some binding sites. Combinatorial polyamine treat- 
ments showed polyamines compete for binding 
to common ribosome sites. Using concepts 
from enzyme kinetics and control engineering, a 
mathematical model of the translational controller 
was developed to describe these complex riboso- 
mal responses to combinatorial polyamine effects. 
Each one of a range of model predictions was suc- 
cessfully validated against experimental frameshift 



frequencies measured in S-adenosylmethionine- 
decarboxylase and antizyme mutants, as well as in 
the wild-type genetic background. 

INTRODUCTION 

Polyamines are essential, ubiquitous polycations found in 
all eukaryotic and most prokaryotic cells. They are utilised 
in a wide range of core cellular processes such as binding 
and stabilising RNA and DNA, mRNA translation, 
ribosome biogenesis, cell proliferation and programmed 
cell death (1). Polyamine depletion results in growth 
arrest (2), whereas their over-abundance is cytotoxic 
(3,4). Thus, homeo statically regulating polyamine 
content within a non-toxic range is a significant regulatory 
challenge for the cell. 

Ornithine decarboxylase (ODC) is the first and rate 
limiting enzyme in the biosynthesis of the polyamines 
(Saccharomyces cerevisiae pathway shown in Figure 1). 
The key regulator of ODC in a wide range of eukaryotes 
is the protein antizyme (5). There is a single antizyme 
isoform in S. cerevisiae, Oazl (6). Antizyme binds to 
and inhibits ODC, and targets it for ubiquitin- 
independent proteolysis by the 26S proteasome (7,8). 
Antizyme synthesis is in turn dependent upon a 
polyamine-stimulated +1 ribosomal frameshift event 
during translation of its mRNA. Antizyme degradation 
by the ubiquitin pathway is inhibited by polyamines (6). 
Polyamines thus regulate their homeostasis via a negative 
feedback system, by controlling Oazl synthesis and 
inhibiting its proteolysis. 

To understand how this negative feedback contributes 
to the complexity of polyamine homeostasis will require a 
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Figure 1. The polyamine biosynthetic pathway in S. cerevisiae. 

(A) Antizyme (Oazl), the main regulator is indicated; it targets ODC 
(Spel) for degradation by the 26S proteasome. Metabolites are shown 
in italics and proteins in Roman type. Asterisks denote genes deleted in 
the spel spe2 paal fmsl deletant strain. See text for further details. 

(B) Block diagram representing the regulation of polyamine biosynthe- 
sis by antizyme as a modular feedback control system. 

robust quantitative characterisation of the underlying 
control mechanisms (9). One partial mathematical model 
of polyamine metabolism has been developed previously 
(10). However, this model lacks a description of antizyme 
regulation through frameshifting, or of polyamine regula- 
tion of antizyme degradation, and is therefore of limited 
use in understanding the consequences of perturbing a 
pathway structured with complex homeostatic controls. 
In control engineering terms, polyamine regulation 
consists of a process (the biosynthetic pathway), which is 
affected by an actuator (the protein antizyme) under the 
control of multiple negative feedback loops (Figure IB). 
In this representation, the translational frameshift event 
plays the role of a feedback controller, and it is clear that a 
complete system model requires a quantitative character- 
isation of this controller. 

In S. cerevisiae, polyamines stimulate antizyme frame- 
shifting at an in-frame UGA codon that terminates 
antizyme ORF1 (6). Frameshifting is enhanced by the im- 
mediately preceding GCG codon, for which there is no 
tRNA whose anticodon matches the codon by Watson- 
Crick or Wobble base pairing in yeast (11,12), and by 
the poor context of the UGA stop codon which causes 
inefficient termination (13,14). Ribosomal frameshifting 
thus occurs in direct competition with an inefficient trans- 
lation termination event and potentially, stop codon 
readthrough by a cellular tRNA. Indeed, polyamines 
have been reported to stimulate stop codon readthrough 
in mammalian systems and Escherichia coli (15,16). Which 
polyamines govern the ratio of frameshifting, termination 
and readthrough at the antizyme recoding site, and the 
concentrations that affect the ratio, is however far from 
clear. Although a number of reports indicate that 



exogenously supplied putrescine, spermidine and 
spermine are all stimulators of antizyme frameshifting, 
these studies were conducted in wild-type cells with 
active polyamine metabolism, making it impossible to 
define which polyamine is the active reagent (17). 

The relatively recent discovery of Oazl in S. cerevisiae, 
and the conserved nature of polyamine metabolism, 
allows this genetically tractable organism to be used as a 
test-bed system to quantify polyamine-induced ribosomal 
frameshifting, and using that analysis, to develop and 
validate a mathematical model of the ribosomal frameshift 
feedback controller. Using a yeast with multiple muta- 
tions, incapable of polyamine synthesis or metabolism, 
the in vivo effects of exogenously supplied polyamines on 
antizyme frameshifting, readthrough and termination 
could be quantified. These measurements were in turn 
used to generate a mathematical model describing 
polyamine effects on ribosome fate at the OAZ1 transla- 
tional recoding site. Successful validation of this model in 
different mutant and gene overexpression backgrounds 
marks a breakthrough in the process of developing a 
complete predictive model of the ribosome as a polyamine 
feedback controller. 



MATERIALS AND METHODS 

Microbial strains and growth conditions 

Yeast strains (Table 1) were grown under standard condi- 
tions on synthetic defined minimal medium [SD; 0.67% w/v 
yeast defined minimal medium without amino acids 
(Formedium), 2% w/v glucose] supplemented with the ap- 
propriate amino acids or nucleotides at 2mg/ml (except 
for leucine, at 6mg/ml). Spermidine (10~ 5 M) was added 
to support growth of polyamine auxotrophs on SD 
medium plates. Escherichia coli strain XLl-Blue (recAl 
endAl gyrA96 thi-1 hsdRH supE44 relAl lac [F' proAB 
lacP ZAM15 TnlO (Tet r )]; Stratagene) was used through- 
out for cloning experiments. 

Generation of a spel spe2 paal fmsl multiple deletant 
yeast strain 

A spel spel paal triple deletant strain was first generated 
by sequentially deleting the entire SPE2 and PAA1 coding 
sequences in a spe l:\kanMX4 background (Y05034; 
Table 1), using short-flanking homology PCR-targeting 
of NatMX4 and HphMX4 markers, respectively (18). The 
spel spe2 paal deletant (IS527; Table 1) was then mated 
with the fmsl deletant strain (Y10595; Table 1), forming a 
diploid that was then sporulated using pre-sporulation and 
sporulation media (19) containing 50 uM spermidine and 
100 uM (3-alanine. Tetrad analysis identified multiply 
deleted spel spe2 paal fmsl strains (IS532; Table 1). 

Determination of in vivo polyamine concentrations 

Polyamines were assayed essentially using the method 
described with the following modifications (20). All glass- 
ware used in polyamine assays was washed with 
concentrated HC1 and rinsed in distilled water to remove 
adherent polyamines. Colonies grown on SD medium 
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Table 1. Saccharomyces cerevisiae strains used in this study 



Strain 


Genotype 


Source 


BY4741 

spel (Y05034) 

spe2 (Y 11743) 

fmsl (Y10595) 

oazl (Y02776) 

spel spe2 paal (IS527/la) 

spel spe2 paal fmsl (TS532/44d) 


MATa his3-l leu2-0 met 15-0 ura3-0 

MATa his3-l leu2-0 metl5-0 ura3-0 spel::kanMX4 (BY4741 genetic background) 
MATa his3-l leu2-0 lys2-0 ura3-0 spe2::kanMX4 (BY4742 genetic background) 
MATa Ius3-1 leu2-0 iys2-0 ura3-0 fmsT:.kanMX4 (BY4742 genetic background) 
MATa his3-l leu2-0 metl5-0 ura3-0 oaz l:\kanMX4 (BY4741 genetic background) 
MATa his3-l leu2-0 metl5-0 ura3-0 speT.:kanMX4 spe2::natMX4 paal::hphMX4 

(BY4741 genetic background) 
MATa his3-l leu2-0 ura3-0 spel::kanMX4 spe2::natMX4 paaT:.hphMX4 fmsT:.kanMX4 


EUROSCARF 
EUROSCARF 
EUROSCARF 
EUROSCARF 
EUROSCARF 
This study 

This study 



agar plates supplemented with 10 _5 M spermidine were 
inoculated into 3 ml SD overnight cultures, and 
sub-cultured overnight to an OD 600 of 0.65-0.85 in 20 ml 
SD medium supplemented with different concentrations of 
polyamines. Cells were harvested, washed twice in ice-cold 
phosphate buffered saline (PBS) and resuspended in 400 ul 
ice-cold PBS (suspension A). To carry out high-pressure 
liquid chromatography (HPLC) analysis, 200 ul of suspen- 
sion A was centrifuged (15 000g, 30 s) and the cell pellet 
resuspended for lh at 4°C in 200 ul of 1.5 M perchloric 
acid, resulting in a volume of 215 ul. Cells (10 ul) were 
enumerated. Cells were then harvested (12 000g, lOmin, 
4°C) and 50 ul of the supernatant diluted in sterile distilled 
water to a final volume of 200 ul. In parallel, polyamine 
standards (Sigma Aldrich) were prepared using putrescine, 
spermidine and spermine over a range of 0.125^1 nmol in 
200 (il of 0.2 M perchloric acid. The commercial prepar- 
ation of putrescine and spermine used contained low-level 
spermidine contamination (<3%), although model fitting 
and data interpretation in every case took account of this 
low-level contamination. Samples were subjected to 
dansylation by alkalinizing with 50 ul of 1 g/ml sodium 
carbonate decahydrate, followed by the addition of 
200 (il of dansyl chloride (15mg/ml in acetone). The 
samples were incubated overnight in the dark at 30° C 
with shaking. Then 125 ul of lOmg/ml proline was 
added and samples incubated at 37°C for 15min. 
Polyamines were extracted using 500 ul of toluene by vor- 
texing for 30 s, and centrifugation (15 000g, 3min). The 
organic phase was recovered and evaporated to dryness 
under nitrogen. Finally, the samples were reconstituted in 
200 ul of methanol and 10 ul of sample subjected to HPLC 
analysis using a Gilson 234 Injector to load samples 
onto a Hichrom HIRBP-2922 column (C8/C18, 5 urn; 
150 x 4.6 mm). Samples were chromatographed using 
a methanol gradient (72-90%; 2min held at 72%, 
increased to 90% >14min, held for 4min, before return- 
ing to 72% and re-equilibrated for 5 min) with a flow of 
1 ml/min. At least three independent transformants of any 
given type were quantified by HPLC. Polyamine intracel- 
lular concentrations assumed a typical yeast cell volume of 
40 urn 3 (21); careful measurements showed the quadruple 
yeast mutant size did not differ from that of the wild type, 
and that addition of polyamines to the growth medium 
left cell size unaffected. 

Plasmid vectors 

To assay ribosomal frameshifting and readthrough at 
the antizyme frameshift site, pairs of complimentary 



oligonucleotides corresponding to the core OAZ1 frame- 
shift sequence plus 17 and 10 nt of the flanking 5' and 3' 
sequences, respectively, were cloned into the unique NotI 
site of vector pAC98T using standard methods (22) 
(Supplementary Figure SI). pAC98T directs the expres- 
sion of a P-galactosidase-luciferase fusion protein and is 
identical to pAC98 vector (23), differing only in the 
promoter that drives expression of the reporter fusion 

(24) . Cloning the pair-wise combination of oligonucleo- 
tides, 5'-GGCCG GGA TTT AAG GAT TGG TGC 
GCG T GAC ATC CCT CTA G-3', 5'-GGCC G GGA 
TTT AAG GAT TGG TGC GCG TGA CAT CCC TCT 
A-3' and 5'-GGCCG GGA TTT AAG GAT TGG TGC 
GCG TGG CAT CCC TCT A-3', into pAC98T produced 
plasmids: pAC98T-OAZl-FS (measuring frameshift 
frequency), pAC98T-OAZl-RT (measuring readthrough 
frequency) and pAC98T-OAZl-Cont, respectively. In 
pAC98T-OAZl-FS, the OAZ1 frameshift site is cloned 
such that luciferase is in the +1 frame relative to lacZ. 
By placing the luciferase ORF in the 0 frame 
(pAC98T-OAZl-RT), downstream of the antizyme 
UGA stop codon, luciferase activity reports stop codon 
readthrough. The remaining plasmid, 
pAC98T-OAZl-Cont, consists of the luciferase ORF in 
the 0 frame, preceded by the OAZ1 frameshift context, 
in which the UGA stop codon was replaced by a TGG 
tryptophan codon. This construct reports 100% produc- 
tion of the P-galactosidase-luciferase fusion. Plasmids 
were transformed into yeast using standard methods 

(25) . High-intracellular concentrations of putrescine were 
engineered using a plasmid expressing a GAL1 
promoter-controlled SPE1 gene (Open Biosystems). SD 
culture media containing a range of galactose concentra- 
tions controlled the expression of SPE1 in pGALl-SPEl 
transformants. Polyamine effects on stop codon 
readthrough in different contexts were assayed as 
described (Supplementary Figure S2, legend) using 
variants of pAC98 in which the lacZ and luc reading 
frames are separated by stop codons in different 
contexts (26). 

Ribosomal frameshifting and readthrough assay methods 

Dicistronic assays for frameshifting and stop codon 
readthrough were performed essentially as described (24) 
with the following modifications. Yeast transformants in 
suspension A (200 ul; see 'Determination of in vivo 
polyamine concentrations' section above) were harvested 
and resuspended in 200 ul lysis buffer (0.1 M potassium 
phosphate buffer pH 7.8, 0.2% Triton X-100) and 
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broken using glass beads (5 mm) in a FastPrep machine 
(MP Biomedicals) over three cycles of 20 s with cooling. 
The lysate was centrifuged (15 000g, lOmin, 4°C) and 
the supernatant diluted 250-fold before assay for 
(3-galactosidase and firefly luciferase in 96 well plates 
using chemiluminescent reporter gene assays, respectively 
Galacto-Light Plus System (Applied Biosystems) and the 
Bright-Glo Luciferase Assay System (Promega) and a 
Glomax 96 luminometer (Promega; 1 s measurement 
after a 2-s delay following injection of 100-ul assay 
reagents). Frameshift frequencies were determined by 
first normalizing measured pAC98T-OAZl-FS luciferase 
activities using the corresponding (3-galactosidase activity. 
This value was expressed as a percentage of the corres- 
ponding normalized luciferase activity measured using 
pAC98T-OAZl-Cont. 

Development of a mathematical model and curve fitting 
to experimental data 

The eight kinetic parameters used in the eventual 
frameshifting function (equation 5) were derived using 
the curve fitting procedure implemented in NLREG 
software (Nonlinear Regression Analysis Program 
Version 6.5, Phillip H. Sherrod, 1991-2010). NLREG 
uses a model/trust-region technique along with an 
adaptive choice of the model Hessian. The algorithm is 
a combination of Gauss-Newton and Levenberg- 
Marquardt methods. 

The foundation of the minimization technique used by 
NLREG is to compute the sum of the squared residuals 
(SSR) for one set of parameter values and then slightly 
change each parameter value and recalculate the SSR to 
analyse how it is affected by the parameter value change. 
By dividing the difference between the initial and updated 
SSR values by the amount the parameter was changed, 
NLREG is able to provide the approximate partial 
derivative with respect to the parameter. This partial 
derivative is employed by NLREG to decide how to 
change the value of the parameter for the next iteration. 

If the initial guess for the parameter is not far from the 
optimum value, the procedure will provide the best 
estimate for the parameter. Such a procedure was 
carried out simultaneously for all parameters representing 
a minimization problem in eight-dimensional space as 
there are eight kinetic parameters in the final frameshifting 
function (equation 5). In this way, the algorithm was used 
to estimate the kinetic parameters that deliver the best fit 
for the frameshifting function and experimental 
measurements. As an important characteristic of non 
linear regression analysis we specifically considered the 
SSR between the experimental frameshifting (y ; ; 
equation 1) and predicted frameshifting (y ; p ; equation 2) 
values at each measured polyamine concentration. 



SSR = J> ; -jf) 2 



(1) 



kinetic parameters. These standard errors can be 
associated with the variance in the experimental dataset 
and allowed us to compute the uncertainty range 
associated with prediction by the frameshifting function 
(equation 5). This uncertainty range is non-linear with 
respect to the magnitudes of polyamine kinetic 
parameters. 

To perform optimization of the polyamine kinetic 
parameters in the frameshifting function (equation 5) by 
means of non linear regression analysis (NLREG 
software) we divided all available data into training and 
testing datasets, the latter to be used in model validation. 
The training data set consisted of a series of single, double 
and triple polyamine effects on frameshift levels in the 
quadruple yeast mutant spel spe2 paal fmsl. To begin 
with we trained the function to predict frameshifting 
from single polyamine data sets. Optimized kinetic 
parameters from the single effect served as an initial 
guess for NLREG analysis of double polyamine effects 
on frameshifting. Two experimental training data sets 
were used for the double polyamine effect on 
frameshifting: the combined spermine and spermidine 
(with no putrescine) data set; the combined spermidine 
and putrescine (with no spermine) data set. The optimised 
parameters from both these data sets are used as two 
initial guesses for the triple polyamine effect on 
frameshifting. At the end stage, both these initial guesses 
deliver a unique optimized polyamine kinetic parameters 
solution expressed by the best fit of the frameshifting 
function (equation 5) to all available frameshift and 
polyamine concentrations in the quadruple deletant 
yeast. The testing data set (the model validation data 
set) consisted of frameshift in oazl, wild-type (WT) 
BY4741, spe2[pGALl-SPEl] and spe2[pRS426] yeast 
strains. We used the frameshifting function with final 
optimized kinetic parameters to predict frameshifting 
levels in the testing data set. 

The single experimental effect of spermidine on 
frameshifting was analysed by non-linear regression 
analysis (NLREG) with the Michaelis-Menten function 
(equation 2) which is a simplified form of the frameshift 
function (equation 5) with zero spermine and putrescine, 



FS([Spd], [Spm] = 0, [Put] = 0) = y™?rc P *?j +*FS 

-Km d +[Spd] 



(2) 



where ^T md designates the Michaelis constant for 
spermidine and K maxd the maximal velocity (% frameshift) 
for spermidine. Spermidine concentration is designated by 
[Spd]. B FS indicates the basal level of frameshifting 
frequency, added as a constant. The best curve fit 
(equation 2) to experimental data has SSR = 5.89 and 
holds the optimal set of kinetic parameters estimated for 
single spermidine effect on frameshifting. 

A similar analysis was performed for the single 
putrescine effect on frameshifting with NLREG model 
of function (equation 3), which is a simplification of 
equation 5 in the absence of spermine and spermidine. 



In addition we also emphasize the regression 
characteristic that is the standard error of the estimated 



V m ax [Putf 

FS([Put], [Spm] = 0, [Spd] = 0) = max ' L v +Bfs 



(3) 
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Here, i^o.sp, designates the affinity constant for putrescine 
and V maX p, the maximal velocity (% frameshift) for 
putrescine. The putrescine concentration is designated by 
[Put]. The Hill coefficient is designated p. B FS indicates the 
basal level frameshifting frequency, added as a constant. 
The best curve fit (equation 3) to experimental data holds 
sum of squared residuals magnitude 44.55 and 
corresponds to the optimal set of kinetic parameters 
estimated for a single putrescine effect on frameshifting. 
The Hill function (equation 4) represents the single 
spermine effect on frameshifting, and its kinetic 
parameters were optimized using NLREG, producing a 
best-fit curve for the single spermine treatment with sum 
of squared residuals of magnitude 2.89. 

V rSnml'" 

FS([Spm], [Spd] = 0, [Put] = 0) = ™ K iJ im +g FS 

(4) 

Here, K Q 5m , designates the affinity constant for spermine 
and Fmaxm, the maximal velocity (% frameshift) for 
spermine. The spermine concentration is designated by 
[Spin]. The Hill coefficient is designated m. B FS indicates 
the basal level frameshifting frequency, added as a 
constant. 

We then analysed the frameshift effect of pairwise 
polyamine treatments on the quadruple yeast deletant to 
further optimise values of the frameshift function kinetic 
parameters. For the pairwise experimental treatments, the 
third polyamine was excluded from the medium. Using 
such data, a search was performed for the best fit of the 
function (equation 5) to predict experimentally observed 
frameshift frequencies. By these means we obtain the 
optimised polyamine kinetic parameters to use as the 
initial guess for triple polyamine effect on frameshifting 
training data set. The NLREG final results for sum of 
squared residues for double spermine and spermidine 
(with no added putrescine) treatment was SSR = 78.18 
and for spermidine and putrescine (with no added 
spermine) was SSR = 118.49. 

RESULTS 

A spel spe2 paal fmsl deletant yeast strain disabled 
for polyamine metabolism 

Addition of any polyamine to the growth medium stimu- 
lates in vivo antizyme frameshifting in cell cultures (17). 
However, these studies were conducted in systems where 
metabolic interconversion of the added polyamines was 
possible, preventing the rigorous analysis of the effects 
of individual polyamines on antizyme frameshifting and 
readthrough. It was therefore essential to define how each 
polyamine governs frameshift and readthrough efficiency 
in a system in which polyamine conversion is prevented. 

To generate such a metabolically isolated system, the 
SPE1, SPE2, PAA1 and FMS1 genes were deleted to 
generate a mutant yeast strain in which de novo synthesis 
and in vivo metabolic interconversion of polyamines taken 
up from the medium is prevented (Figure 1). The lack of 
Spel and Spe2 prevents putrescine synthesis and its 



conversion to spermidine and spermine (27,28,29,30), 
while FMS1 deletion prevents the reverse spermine to 
spermidine reaction (31,32,33). The absence of Paal 
prevents polyamines from being acetylated (34) (Figure 
1). Polyamine analysis of the now polyamine auxotrophic 
quadruple deletant confirmed that de novo synthesis had 
been eliminated (Table 2). A block in metabolic inter- 
conversion of polyamines was confirmed by separately sup- 
plying each of the three polyamines exogenously to the 
deletant (10~ 2 mM), and subsequently detecting only 
that single polyamine in cell extracts (Supplementary 
Table SI). 

Putrescine, spermidine and spermine stimulate antizyme 
frameshifting in qualitatively and quantitatively 
different ways 

In order to understand how polyamines govern yeast 
antizyme synthesis, we analysed in vivo how each of the 
different polyamines affect the three mutually exclusive 
processes of frameshifting, readthrough and termination 
at the antizyme frameshift site in the quadruple deletant 
yeast. Frameshift and readthrough were measured in 
actively growing wild-type and spel spe2 paal fmsl 
cultures using dicistronic reporter plasmids containing 
the antizyme frameshift site, along with readthrough 
controls ('Materials and Methods' section). In the 
absence of added polyamines, ~16.1% of the translating 
ribosomes shifted to the + 1 frame in the wild-type genetic 
background, and only 3.7% in the quadruple deletant 
(Figure 2). The readthrough frequencies for both 
wild-type and quadruple deletant strains were similar 
and low (0.52 and 0.43%, respectively; Figure 2). Thus 
frameshifting, but not stop codon readthrough, appears 
to be polyamine sensitive. 

The effect of polyamines on frameshifting at the OAZ1 
frameshift site was next examined in the quadruple 
deletant strain grown in media with varying extracellular 
polyamine concentrations. The addition of increasing 
extracellular concentrations of putrescine (10 -1 , 1, and 
lOmM) produced either no effect, or a very minor 
1.7-fold stimulatory effect on frameshift frequencies 
(Figure 3A). However, other studies have shown that 
very high concentrations of putrescine can stimulate 
frameshifting at the yeast Tyl retrotransposon frameshift 
site, specifically under conditions where very low concen- 
trations of spermidine and spermine were present (35). We 
therefore overexpressed the SPE1 gene under the control 



Table 2. Intracellular polyamine concentrations wild-type BY4741 
and spel spel paal fmsl deletant strains' 1 



Polyamine 




Strains 


BY4741 


spel spel paal fmsl 


Putrescine (mM) 


0.17 ± 0.01 


ND b 


Spermidine (mM) 


2.88 ± 0.53 


ND b 


Spermine (mM) 


0.74 ± 0.15 


ND b 



"Analyzed by high-pressure liquid chromatography (HPLC). 
b ND, not detectable by HPLC (detection limit was 0.02 mM). 
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Wild-type spel spe2 paal fmsl 

Figure 2. Antizyme ribosomal frameshifting and readthrough in 
wild-type BY4741 and spel spel paal fmsl deletant strains. Average 
percent frameshifting (FS) and readthrough (RT) for wild-type and 
quadruple deletant strains grown in polyamine-free media. Filled bars 
represent average frameshift and open bars represent readthrough. 
Standard deviations are presented for all bars (typically <15% of 
mean value; n = 3). 



of the GAL1 promoter to engineer higher intracellular 
concentrations of putrescine. Using galactose in the 
culture medium to induce SPE1, extremely high levels of 
putrescine (18.6-71 mM) were attained. At these very 
high-intracellular concentrations, frameshifting was 
indeed stimulated, to a maximum of 34% (Figure 3B). 
Strikingly, the putrescine frameshift response curve was 
sigmoidal, indicating that the binding of putrescine to 
the ribosome to stimulate frameshifting was in some way 
cooperative and may involve multiple binding sites 
(Figure 3B). 

Treating the quadruple deletant cultures with a range of 
extracellular spermidine concentrations (5 x 10~ 4 , 10~ 3 , 
2.5 x 10~ 3 , 10~ 2 , 10 -1 and 1 mM) produced a progressive 
increase in frameshifting, reaching a maximum of 22.0% 
with an intracellular concentration of spermidine of 
8.36 mM (Figure 3C). Spermidine is thus a strong stimu- 
lator of antizyme frameshifting. Finally, treatment of the 
deletant with increasing extracellular concentrations of 
spermine (10"\ 10" 2 , 10" 1 , 0.8, 0.9 and 1 mM) produced 
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Figure 3. Single treatment effects of putrescine, spermidine and spermine on frameshifting at the antizyme frameshift site. Average percent 
frameshifting (filled circles) was measured using a dicistronic assay and plotted versus intracellular polyamine intracellular concentrations in the 
spel spe2 paal fmsl deletant strain. (A) Putrescine effects on frameshifting. The highest putrescine concentration used contained 0.16mM 
contaminating spermidine, and was therefore not used in the curve fitting process. (B) Putrescine-stimulated frameshift frequencies were measured 
in the quadruple deletant strain transformed with pGALl-SPEl grown on a range of galactose concentrations to regulate SPE1 expression. (C) The 
effect of intracellular spermidine on frameshift frequency. (D) The effect of intracellular spermine on frameshift frequency. The highest concentration 
contained 0.2 mM contaminating spermidine. For reference, the wild-type strain BY4741 frameshift frequency (open circles) is represented on all 
graphs. Error bars (horizontal and vertical) indicate standard deviations for three independent transformants, analysed in triplicate. 
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similar increases in frameshifting. A maximal frameshift 
stimulation of 27.6% (a 7.5-fold increase) was observed at 
an intracellular concentration of 2.67 mM. Spermine is 
thus an even more potent stimulator of frameshifting 
than spermidine. The sigmoidal nature of the spermine 
response curve also suggests that binding of this 
polyamine to the ribosome is co-operative and/or 
involves multiple binding sites (Figure 3D). 

Polyamines do not stimulate stop codon readthrough 
in yeast 

In contrast to their marked effect on antizyme frameshift- 
ing, polyamines trigger only slight increases (1.5-fold) in 
readthrough of the antizyme premature UGA stop codon 
in yeast (Supplementary Figure S2A, B and C). These 
findings contrast with other studies reporting that 
polyamines can enhance stop codon readthrough in 
E. coli and in mammalian systems (15,16), and enhance 
readthrough between 2- and 9-fold at the antizyme 2 
premature stop codon in COS7 cells (17). These observa- 
tions call into question whether polyamines in yeast affect 
frameshifting specifically, or whether they more generally 
enhance mis-translation events. Accordingly, the effect 
of polyamines on translational readthrough of stop 
codons in a variety of different contexts was assessed 
using dicistronic reporters in the quadruple deletant 
yeast strain. Exogenous polyamine concentrations were 
chosen that significantly enhance translational frameshift 
frequencies (see legend, Supplementary Figure S3). The 
results show clearly that irrespective of the nature of the 
stop codon and its context, neither spermidine nor 
spermine significantly enhanced readthrough of any of 
the three stop codons (Supplementary Figure S3). 

Frameshift stimulation by polyamine combinations reveals 
evidence of competition for common ribosomal binding 

sites 

Having established that spermidine and spermine (and at 
high concentrations, putrescine), stimulate antizyme 
frameshifting, combination-specific polyamine effects on 
frameshifting and readthrough were tested. The spel 
spe2 paal fmsl deletant strain was therefore cultured in 
medium containing different combinations of extracellular 
polyamines (spermidine and spermine; putrescine and 
spermidine; putrescine and spermine). The results 
showed that together, spermidine and spermine 
(combined 0.01-1.04mM) produced approximately 
median frameshift frequencies between the values that 
each polyamine stimulated individually. At higher 
combined concentrations of spermidine and spermine 
(>2.5mM) frameshift frequencies were observed higher 
than those obtained when testing spermidine and 
spermine individually (Figure 4A). 

The effect of putrescine, which at low concentrations 
does not stimulate frameshifting (Figure 3A), was then 
assessed in combination with either spermidine or 
spermine, which individually were strong frameshift 
stimulators. Across a range of putrescine concentrations, 
the stimulatory effects of adding either 0.0025, 0.01 or 
1 mM spermidine were somewhat suppressed in 




Figure 4. Combinatorial effects of polyamines on frameshifting at the 
antizyme frameshift site; experimental data and model prediction. 
Frameshifting frequencies were measured using a dicistronic reporter 
vector containing the OAZl recoding site. Intracellular polyamine con- 
centrations were measured by HPLC. The combined effect of spermi- 
dine and spermine (A), putrescine and spermidine (B) and putrescine 
and spermine (C) on frameshifting are plotted as filled circles, and the 
model prediction (mesh surface) depicting the fitted frameshift function 
(Equation 5). Note that panel (C) data points are derived from 
combined putrescine/spermine treatments some of which contained 
trace spermidine contamination, but the mesh surface indicates pre- 
dicted responses to dual spermine/putrescine treatment only. For 
clarity, experimental data standard deviations are not presented but 
were typically <15% of mean value (« = 3). 
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comparison to the expected levels of frameshifting 
generated by spermidine alone (Figure 4B). However, 
when the level of putrescine was maintained, but the con- 
centration of spermidine was incrementally raised, this 
minor inhibitory effect of putrescine was overcome 
(Figure 4B). No such suppression of spermine-induced 
frameshifting by putrescine was observed (Figure 4C). 
At low concentrations, putrescine therefore acted as a 
weak suppressor of the stimulatory effect of spermidine 
on frameshifting. Polyamines thus do not significantly 
stimulate stop codon readthrough, but exert combinator- 
ial and competitive effects on frameshifting consistent 
with their similar chemical nature. 

Modeling the l iposomal frameshift feedback controller 

Using the experimental data above, a mathematical 
function was introduced to capture the complex individual 
and combinatorial effects of the three polyamines on ribo- 
somal frame maintenance. We noted that the single treat- 
ment spermidine frameshift response curve (Figure 3C) 
bore a striking similarity to a Michaelis-Menten enzyme 
kinetic function. Frameshifting was thus represented as an 
enzyme-like activity of the ribosome, stimulated by 
spermidine 'substrate'. After first subtracting the back- 
ground frameshift frequency from all measurements, a 
curve-fitting analysis was conducted to determine the 
maximal frameshifting (V max ), and affinity constant 
(K m ) for spermidine interaction with the ribosome. 
Frameshift response curves for single putrescine and 
spermine treatments had indicated co-operative binding 
effects (Figures 3B and D), so Hill functions were used 
to represent their individual effects on frameshifting, and 
curve fitting (using non-linear regression methods; 
Materials and Methods) used to minimize residual sum 
of squares and define the optimal frameshifting (K max ), 
affinity constant (K Q5 ), and Hill coefficients for 
ribosome interaction (Table 3). 

To describe the combinatorial effects of the three 
polyamines, we recognized the chemical and structural 
similarities of the polyamines and proposed that they act 
as competing pseudo-substrates for the same active site(s) 
on the ribosome. Consistent with this, spermine and 
spermidine combinations can produce intermediate frame- 
shift frequencies at certain concentrations (Figure 4), con- 
sistent with their structural similarity. We therefore 
generated an enzyme kinetic function 'FS' based on the 
kinetics of an enzyme with competing substrates (36) that 



described the action of the pseudo-substrates (Put, Spd, 
Spm) on a single enzyme activity (i.e. frameshift effi- 
ciency). The total frameshift frequency was, therefore, 
the sum of three individual frequencies contributed by 
each of the three polyamines. In the global frameshift 
function (equation 5), subsidiary frameshift functions 
describe the contribution made by each polyamine to the 
total frameshift frequency, moderated by competition 
from the other two polyamines through additional [con- 
centration]/[affmity constant] terms in the denominator of 
each subsidiary term as described for an enzyme with 
competing substrates (36); 

FS([Put],[Spd],[Spm]) 

KmaxJPut]" 



K maXd [Spd] 



^mdl 1 



. [Spm]"' , [Put]" 

K m K p 

[Spm] 



+[Spd] 



(5) 



' max B 



FS 



In equation 5, 'K' parameters define affinity 
constants, 'V parameters indicate velocities (frameshift- 
ing) (see Table 3 for details) and concentrations of the 
individual polyamines are indicated by [Spd], [Spm] and 
[Put]. The Hill coefficients are designated by m (spermine) 
and p (putrescine). B FS indicates the constant basal-level 
FS measured in the absence of polyamines. 

The feedback-controller function describing the com- 
peting polyamine substrates accurately predicts the 
combinatorial effects of pair-wise polyamine treatments, 
which accurately map onto the 3D model mesh (Figure 4), 
validating our assumption that the three polyamines ef- 
fectively compete for the same site(s) on the ribosomal 
'enzyme'. The biological principle underpinning the math- 
ematical model is thus an accurate representation of the 
polyamine molecular interactions with the ribosome. 

Validation of the frameshift feedback controller model 
using polyamine pathway mutants 

To rigorously validate the model of the ribosomal frame- 
shift feedback controller, other engineered yeast genetic 
backgrounds were used to generate hitherto unexplored 



Table 3. Polyamine frameshift model kinetic parameters 



Polyamine 


Kinetic parameter; symbol 




Value 


Standard error 
(single polyamine) 


Uncertainty range 
(combined treatment) 


Putrescine 


Maximum velocity K maX(i (percentage 


of frameshift) 


32% 


±7.6 


24.4-39.6 




Affinity constant: A^ 5p 




26.6 mM 


±7.4 


19.2-33.9 




Hill coefficient: 




2.4 


±1.4 


1.1-3.8 


Spermidine 


Maximum velocity K maXd (percentage 


of frameshift) 


33.9% 


± 12.5 


21.4-46.3 




Micliaelis constant: K md 




6.8 mM 


±4.2 


2.6-10.9 


Spermine 


Maximum velocity K maXm (percentage 


of frameshift) 


21.3% 


±2.7 


18.6-24 




Affinity constant: K 05nl 




1.7 mM 


±0.2 


1.5-1.8 




Hill coefficient: m 




5.2 


±2.3 


2.9-7.5 
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combinations and concentrations of the three polyamines. 
Accordingly, frameshifting and intracellular polyamine 
concentrations were measured in a spe2 mutant with and 
without SPE1 gene overexpression, in an oazl deletant, 
and in the wild-type strain. The measured intracellular 
polyamine concentrations from each mutant strain were 
fed into the frameshift function, and the predicted in 
silico frameshift values compared with the corresponding 
in vivo measurements. 

In the spe2 deletant strain, putrescine accumulated 
to moderately high levels (6.5 mM; Figure 5A). The 
measured 4.5% frameshifting was accurately predicted 
by the model (4.7% frameshift, Figure 5B). Over- 
expressing the SPE1 gene under the control of the GAL1 
promoter in this strain (spe2 pSPEl), produced 65.3 mM 
intracellular putrescine and a frameshift frequency of 
27.4%, again, well approximated by the model (32.4%; 
Figure 5B). To further validate the model, the oazl 
antizyme deletant was selected, which generates atypical 
concentrations of all three polyamines (high levels of 
putrescine and spermidine, normal spermine; Figure 5A), 
and high levels of frameshifting (21.7%), and which thus 
tested the ability of the model to integrate the competing 
polyamine effects. This level of frameshift was accurately 
predicted by the feedback controller function (21.4%; 
Figure 5). Finally, we verified that the wild-type frameshift 
frequency (16.1%) was accurately predicted by the model 
(13.9%; Figure 5B). The model thus demonstrated an 
ability to accurately predict the frameshift efficiencies 
generated by a wide range of different genetic conditions 
that were typified by markedly different polyamine steady 
states. 



DISCUSSION 

In this study, a novel quadruple yeast gene knockout 
strategy was devised in which de novo synthesis and meta- 
bolic interconversion of supplied polyamines is prevented. 
Using this tool, this study aimed to define for the first 
time how each of the polyamines individually stimulates 
the frameshifting, termination and readthrough events at 
the antizyme frameshift site in vivo. With this information, 
a predictive mathematical model of the ribosomal frame- 
shift controller (Figure IB) could then be generated which 
will form the cornerstone of future systems biological 
analysis of polyamine biosynthesis. 

The observed polyamine effects on antizyme stop codon 
readthrough in the quadruple knockout yeast strain 
(Supplementary Figure S2) showed clear, system-specific 
differences in the way that polyamines influence translation 
termination. Previous studies have reported stimulation of 
stop codon readthrough in E. coli (16). Polyamines can also 
stimulate readthrough in in vitro translation systems (15). 
However, in this study, readthrough of the antizyme stop 
codon, and of other stop codons in a range of weak and 
strong contexts was refractory to polyamine stimulation 
(Supplementary Figures S2 and S3). This highlights poten- 
tial differences between yeast and E. coli ribosome re- 
sponses, and between in vivo measurements and those 
made in cell-free systems. However, the observations 




spe2 spe2 pSPEl oazl Wild-type 




spe2 spe2 pSPEl oazl Wild-type 



Figure 5. Validation of the ribosomal controller frameshift model. 
(A) Polyamine concentrations were measured in a spe2 mutant with 
and without SPEl gene overexpression, in an oazl deletant, and in 
the wild-type strain (putrescine; filled bar, spermidine; light grey, 
spermine; dark grey). Error bars represent standard deviations 
(n - 3). (B) Ribosomal frameshift frequencies were measured in the 
same mutant panel (filled bars) and the intracellular polyamine concen- 
trations measured in the mutants were fed into the frameshift function 
to predict the frameshift frequency (open bars; error bars represent 
model uncertainty range originating with variation in the experimental 
data). 



made simplify the development of a mathematical descrip- 
tion of the yeast ribosomal feedback controller, since only a 
two-way competition between frameshifting and termin- 
ation need be considered. 

The analysis of polyamine frameshift responses revealed 
that although spermidine stimulates frameshifting with 
a hyperbolic (Michaelis-Menten) function, in contrast, 
putrescine and spermine each appear to bind to the 
ribosome in a co-operative manner (Figure 3B and D; 
Hill coefficients of 2.4 and 5.2, respectively). We speculate 
that this co-operativity could be generated by the existence 
of multiple polyamine binding sites on the ribosome, and 
that frameshifting stimulation requires gradual polyamine 
occupancy of these sites. This conclusion is consistent with 
photo-affinity labeling studies revealing spermine binding 
at several sites in the 23S rRNA of E. coli, including 
those regions involved in positioning and translocation 
of tRNA (37,38,39). Why spermidine responses are not 
co-operative is unclear (Figure 5B), but by extension, 
this may suggest the existence of a single spermidine 
binding site. This could in turn explain why at low con- 
centrations putrescine slightly inhibits spermidine-induced 
frameshifting; putrescine occupancy of the putative, single 
spermidine-sensitive site may compete with spermidine 
binding, without putrescine itself triggering frameshifting 
(which requires multiple site occupancy; Figure 3B). 
However, the high-affinity constant for putrescine- 
ribosome interaction (10-fold higher than spermidine; 
Table 3) ensures that this is only a subtle effect. 
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Putrescine stimulation of antizyme frameshifting has 
been reported previously, although in cell culture 
systems where metabolic interconversion of putrescine 
into spermidine and spermine could occur (17,40). In 
our study, where interconversion was prevented by gene 
knockout, low concentrations of putrescine did not stimu- 
late frameshifting (Figure 3A). This would be consistent 
with the fact that putrescine is a pathway intermediate 
with limited polycation character and function; as such, 
it is not desirable for low concentrations of putrescine to 
trigger frameshifting, antizyme production and pathway 
shut-down. Previously, reports indicated that a yeast 
spe2 deletion (exhibiting somewhat raised putrescine con- 
centrations) triggered high levels of frameshifting at 
the yeast retrotransposon Tyl frameshift site, but only 
following a period of spermidine deprivation (35,41). 
Experimentally, our study also confirmed that putrescine 
stimulates frameshifting at high concentrations 
(Figure 3B). Our experimental data, and the derived math- 
ematical frameshift function, also indicates that spermi- 
dine would suppress the putrescine-induced frameshifting 
to some extent (Figure 4B). It is possible that any differ- 
ences between our study and others (35,41) in the degree 
of spermidine suppression of putrescine frameshifting 
might stem from their use of a Tyl frameshift site, 
rather than antizyme, to report polyamine responses. 

Pair- wise treatments of polyamines (Figure 4) revealed 
that spermidine and spermine together generate frameshift 
frequencies intermediate between their individual effects. 
Consistent with their similar positively charged chemical 
structure, this indicated that polyamines might compete 
for the same binding sites on the ribosome. We used 
this discovery as the platform for the development of a 
mathematical model capable of describing frameshift re- 
sponses to a full range of polyamine conditions in the cell. 
Importantly, although the model parameter fitting process 
primarily made use of single and pair-wise polyamine data 
sets, the complete model was able to accurately predict 
frameshift efficiencies measured in both the wild-type 
strain, and in an antizyme mutant (oazl), each of which 
contain very different (triple) polyamine combinations. 

The successful process of validating model predictions 
supports the underlying assumption that, in vivo, 
polyamines compete for common polyamine binding 
sites on the ribosome, and thus individually act as inhibi- 
tors of frameshift stimulation by either of the other two 
polyamines. The approach has also allowed the in vivo 
parameterisation of that competition. The distinct hyper- 
bolic and co-operative responses of frameshifting to 
spermidine and spermine respectively determine non- 
linear frameshift responses to the polyamine milieu. This 
regulatory complexity justifies the need for a mathematical 
function to describe the four-dimensional relationship 
between frameshifting and the three polyamines. The 
model in addition provides a robust and quantitative 
tool to analyze the action of therapeutic polyamine ana- 
logues and other potential frameshift agonists on transla- 
tional frameshifting. The validated ribosomal frameshift 
model will in addition act as a fulcrum for the develop- 
ment of complete mathematical models of the polyamine 
pathway and its regulation. 



Beyond the field of polyamine biology, the ribosomal 
controller model we have developed is directly relevant to 
understanding other feedback translational regulation 
systems, for example the CPA1 arginine-regulated expres- 
sion system (42). Translational feedback modules, with 
their inherently rapid response characteristics, could also 
represent useful building blocks in synthetic biology. In 
contrast to the well-understood example of transcriptional 
feedback, the current predictive mathematical understand- 
ing of translational control is very poor, representing a 
major stumbling block to its use in the design of synthetic 
circuits. The model presented in this paper therefore not 
only helps define the workings of a key metabolic pathway 
important in a number of human pathologies, but will also 
allow the introduction of novel translational-based 
feedback control modules to the synthetic biologist's 
toolkit. 
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